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Summary 

The subsonic flow field about the basic F-16XL air- 
plane model, with and without air dams, was analyzed by 
investigators using an inviscid unstructured grid tech- 
nique. Computed surface pressure distributions at 13 
fuselage stations and 10 butt-line stations were compared 
with a wind-tunnel investigation at Mach 0.148 for a 
range of angles of attack from 0° to 20°. To evaluate the 
effect of grid dependency on the solution, a grid study 
was performed in which fine, medium, and coarse grid 
meshes were generated for the F-16XL without air dams. 
The off-surface grid field, which was locally adapted to 
the vortical flow field and compared to the nonadapted 
flow field, showed improved correlation with the wind- 
tunnel data. The computed off-body flow quantities for 
the no-air-dam configuration are compared to five-hole 
pressure probe data at a = 10°. The grid for the basic 
F-16XL with air dams also was adapted to the vortical 
flow field, and solutions were obtained at a = 10°, 13°, 
and 15°. A comprehensive analysis of the off-body com- 
puted pressure contours and velocity vectors is presented 
for configurations with and without the air dam. 

Introduction 

This paper describes a flight-wind-tunnel- 
computational fluid dynamics (CFD) validation experi- 
ment for the F-16XL airplane (ref. 1 and work in 
progress at Langley Research Center by S. J. Rickard). 
This paper focuses on the application and calibration of 
an unstructured grid method to solve the Euler equations 
for 0.04-scale models tested in the Basic Aerodynamic 
Research Tunnel (BART) facility (ref. 2) at Langley 
Research Center (fig. 1). The use of unstructured grids to 
compute flows about complex geometries has increased 
in recent years because the unstructured tetrahedral mesh 
discretizes irregularly shaped domains more easily than a 
structured grid does. Given the drastically reduced time 
to generate a grid, this approach is especially suitable for 
use in the design stages of configurations. Two disadvan- 
tages of using unstructured grids encountered by the 
investigator were increased computer memory require- 
ments and increased postprocessing time. Solving the 
Euler equations requires less central processing unit 
(CPU) time than do the more accurate Navier-Stokes 
equations, and depending on the flow characteristics, the 
Euler equations can offer reasonable results fairly 
quickly. 

One obvious shortcoming to using an Euler solver is 
the exclusion of viscous flow effects. While the Navier- 
Stokes equations are necessary to accurately capture all 
the flow physics, researchers have found that many of the 
important features of the primary vortex and its interac- 
tion with the wing can be modeled by the Euler equa- 


tions. Euler equations work well for vortex-dominated 
flows generated by slender delta wings, where the vis- 
cous effects do not dominate the off-surface flow field 
and the flow separates from a sharp leading edge (ref. 3). 
The sensitivity of the Euler codes to such numerical 
parameters as artificial viscosity of the numerical algo- 
rithm and discretization errors also affects the solution 
(ref. 4). Therefore, the applicability of the Euler equa- 
tions to this type of problem will depend on these types 
of considerations. Recent unstructured grid studies per- 
formed with the unstructured grid Euler code USM3D 
on a wing-pylon-store configuration (ref. 5) and an 
isolated fuselage geometry (ref. 6) reveal favorable com- 
parisons of surface pressure distributions as well as 
forces and moments. Because favorable results were 
obtained from previous studies and the F-16XL has a 
slender cranked-arrow wing, USM3D was used to ana- 
lyze this configuration. 

Symbols and Abbreviations 

BL butt-line stations 

CFD computational fluid dynamics 

CFL Courant-Friedrichs-Lewy number 

CPU central processing unit 

C L lift coefficient 

Cp static pressure coefficient 

Cp t total pressure coefficient 

FS fuselage station 

M free- stream Mach number 

MW megaword 

RN Reynolds number 

b wing span, 15.58 in. 

c local chord 

c mean aerodynamic chord, 11.87 in. 

nd without air dam 

wd with air dam 

x/c fraction of local chord 

y span wise distance, in. 

y/(b/2)i fraction of local semispan; note local semi- 
span does not include missile rail and missile 

z normal distance, in. 

a angle of attack, deg 

Model and Test 

A three-view schematic of the 0.04-scale basic 
F-16XL configuration appears in figure 2. The basic 
model features an “S-shaped” leading edge at the wing 



apex as opposed to a straight leading edge (modified 
apex), which was also experimentally tested but not 
included in this investigation. The model is 25.41 in. 
long with a 15.55-in. wingspan. The wing has an inboard 
leading-edge sweep angle of 70° and an outboard 
leading-edge sweep angle of 50°. Other model features 
include a vertical tail and a drag chute housing. An air 
dam on each wing extends part way onto the upper 
surface of the aileron-actuator pod. A missile consisting 
of four front and four rear fins attaches to each wingtip 
via a missile rail. The model is constructed with a flow- 
through inlet and nozzle. 

Investigators obtained the experimental data in the 
BART facility (fig. 1) by testing a 0.04-scale F-16XL 
over a range of test conditions: 0° to 20°, 0.07 < M 
< 0.165, and 500,000 < RN < 1 120 000, all at zero side- 
slip. Data types taken included surface static pressure 
(ports mapped in fig. 3) and oil flow, five-hole pressure 
probe in the flow field (for total pressure and velocity 
magnitude and direction), vapor screen, and pressure- 
sensitive paint on several test configurations. Most of the 
data are reported in reference 1 and in the work in 
progress by S. J. Rickard at Langley. 

Selected data samples for the basic model with mis- 
siles and with and without air dams are used in this paper 
for comparison purposes. Geometrical consistency of 
model and CFD numerical surface description are 
assumed in this data comparison. 

Computational Grids and Method 

The surface and volume grids were constructed with 
the grid generator VGRID (ref. 7), which is a tetrahedral 
unstructured grid generator based on the advancing front 
method (refs. 8 and 9). The parameters for the initial 
front (initial surface grid) were defined in GRIDTOOL 
(ref. 10). The GRIDTOOL program allows the user to 
interactively divide the configuration into bilinear 
patches, as well as define the node (or point) and line 
source locations (fig. 4), magnitudes, and directions. The 
nodal and line sources control the spacing on the surface 
and the cell sizes in the grid field. No grid stretching was 
used; a cell aspect ratio of 1 .0 was specified, although not 
always maintained by the code. An initial front was cre- 
ated and projected onto the original database surface with 
GRIDTOOL. The front was then advanced into the field 
by adding tetrahedral cells to complete the grid. 

Unstructured grids were generated for the basic 
F-16XL model with and without air dams. In an attempt 
to minimize the total number of cells, thereby reducing 
memory and run time required, fine, medium, and coarse 
grids of the F-16XL configuration without air dams were 
generated first. This method permitted not only assess- 
ment of the effects of grid dependency on the solution 


but also determination of the most efficient grid in terms 
of solution accuracy and computer resources. 

Grid Study 

An unstructured fine grid was generated first for the 
basic F-16XL without air dams. To ease the grid con- 
struction process, a small modification was made to the 
four fins located at the rear of the missile; the thickness 
of the leading and trailing edges of the fins was 
increased. Also, unlike the database, which modeled the 
tips of the four rear fins as collapsed lines, the computa- 
tionally defined fins were modified to have a finite thick- 
ness. It would have been advantageous to model the fin 
tips with a singular line, which would result in fewer 
total grid cells being required because of the absence of 
the small tip surface. However, the collapsed fin-tip 
geometry produced difficulties for the advancing-front 
method. 

The surface of the missile, missile rail, and surround- 
ing area required a dense distribution of points to resolve 
the small fin geometry and to capture the flow physics in 
this region. A dense point distribution also had to be 
maintained in this area to ensure a good (not skewed) 
local meshing of the cells. Figure 5(a) shows the surface 
and reflection plane of the F-16XL unstructured fine 
grid. The grid outer boundaries were located about 5 c 
from the surface in all directions. The complete fine grid 
is composed of 1 144077 cells and 205317 points, with 
19000 points on the surface. Figure 5(b) shows a close- 
up view of the surface grid on the missile and missile 
rail. As the figure shows, small cell spacing over a fairly 
large field area was necessary to ensure local and global 
meshing with a minimum amount of skewness. A coarser 
grid density was generated by increasing all the source 
strengths by 1.25. A total of 899083 cells and 161967 
points comprised this grid, with 16351 grid points on the 
surface. Likewise, the coarse grid was generated by 
applying a 1.50 increase in the source strengths for the 
fine grid. The resulting coarse grid consisted of 532622 
cells and 96835 points, with 1 1 740 points on the surface. 

The fine, medium, and coarse surface grids are 
shown and compared to one another in figure 6. While a 
gradual decrease in point density is seen from fine to 
coarse grid, higher grid point density is consistently 
maintained for all three grids on the leading edge, in the 
missile area (including the wingtip), on the actuator pod, 
and at the wing/fuselage interface. The off-surface grid 
densities for each grid appear in figures 7 and 8 at fuse- 
lage stations FS10.0 and FS18.0, respectively. The dif- 
ferent grid densities are not as apparent at FS18.0 as they 
are at FS10.0 because of the very small spacing initially 
required to resolve the fins on the fine grid. Even after 
increasing all point and line source strengths by a factor 
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of 1.5, the spacing on and around the missile region 
remained quite small. 

The medium surface grid definition was chosen and 
then locally adapted to the vortical flow regions at 
M - 0.148 and a = 10°. The justification for choosing 
the medium grid was based on solution accuracy and 
computational efficiency and is discussed in detail in the 
“Results” section. The “Results” section also demon- 
strates the need for local grid adaption in the regions of 
vortical flow. In order to find the location of the vortical 
flow areas, the pressure coefficient C P contours from the 
converged solution were plotted in several crossflow 
planes so that the vortical cores (identified by the lowest 
C P value) could be located. Line sources placed through 
these points produced a higher density of grid cells in the 
localized area of the vortex paths. These off-surface line 
sources also affected the point distribution on the surface 
of the adapted medium grid, as shown in figure 9. The 
completed adapted grid consisted of a total of 826420 
cells and 149299 points, with 16132 grid points on the 
surface. The researcher found that by maintaining 
approximately the same number of surface points as the 
medium grid (16351 medium grid surface points), the 
volume grid could be generated more efficiently, cluster- 
ing points only in the vortical regions. As will be shown, 
the adapted grid yielded significantly better surface pres- 
sure correlations with experiment, using fewer grid 
points. Cross-sectional planes at FS10.0 and FS18.0 for 
both the medium and adapted medium grid appear in fig- 
ures 10 and 11, respectively. Note the clustering of cells 
above the wing, which represents the primary vortex 
regions [figs. 10(b) and 11(b)]. The investigator believed 
that the high density of points in the vicinity of the mis- 
sile and wingtip generated by the fine spacing was not 
needed for solution accuracy, and specified a less clus- 
tered field grid around the missile for the adapted 
medium grid (fig. 11). 

Once the investigator gained confidence in unstruc- 
tured grid generation for this configuration, the air dam 
was added to the fine grid using GRIDTOOL. The origi- 
nal database of the air dam had a constant thickness; 
however, to facilitate gridding, the top of the air dam was 
modeled with zero thickness. Additional sources had to 
be placed on the air dam to resolve the geometry and the 
channel flow in this region. The initial front and volume 
were generated by increasing the source strengths by 
1.25, the same source strength as that used for the 
medium grid. The volume grid was adapted for a = 10° 
by using the same source specifications used for the no- 
air-dam configuration at a = 10° because the flow field 
was not expected to change much. The resulting grid 
contained 1052737 cells and 189283 points, with 18392 
grid points on the surface. Figure 12 shows the surface 
grid of this configuration, as well as a close-up view of 


the air-dam geometry and point distributions. As the 
“Results” section shows, the grid adapted for the a = 10° 
solution proved to be sufficient for the a = 13° solution, 
so a new grid did not have to be generated. However, a 
new adapted grid had to be constructed for the a = 15° 
case because the vortical core path had moved far enough 
outside the densely clustered area used for the a = 10° 
through 13° range. A newly generated grid reflected the 
new locations of the vortical core path. The researcher 
determined this new path by locating the vortical core in 
several crossflow planes obtained from a partially con- 
verged solution (to minimize use of computer resources) 
at a = 15° using the grid adapted for the lower angle-of- 
attack solutions. 

Computational Method 

The flow field was computed using USM3D 
(ref. 11), which solves the time-dependent three- 
dimensional Euler equations in a computational domain 
discretized by tetrahedral-mesh elements. Spatial discret- 
ization was accomplished by using the flux-splitting 
method of Roe, which is based on a cell-centered finite- 
volume approach. The solution was advanced in time 
through use of an implicit Gauss-Seidel scheme, and 
convergence was accelerated to steady state by local time 
stepping and implicit residual smoothing. Flow tangency 
was imposed on solid boundary surfaces, and density and 
pressure boundary conditions were set to the cell- 
centered value. Characteristic boundary conditions were 
applied to the far-field subsonic boundary in which the 
fixed or extrapolated Riemann invariants were used, 
depending on the wave direction. Further details about 
the computational method are discussed in reference 1 1 . 

Results 

The solution results for each grid are presented and 
discussed in the same order as outlined in the “Computa- 
tional Grids and Method” section, namely grid study on 
F-16XL model (no air dam), adapted grid for a = 10°, 
and adapted grid with air dams for a = 10°, 13°, and 15°. 
All solutions obtained were at M = 0.148. 

Convergence and Performance Characteristics 

The computations performed in this study (including 
the grid generation) were done on the Cray-C90 located 
at the Ames Numerical Aerodynamic Simulation (NAS) 
center and the Cray-YMP computer located at Langley 
Research Center. Figure 13 shows the convergence histo- 
ries for the coarse, medium, and fine grids obtained at 
a = 10° and M = 0.148. The fine grid required a 
206-MW memory allotment and took 1650 iterations to 
reduce the residuals by 2.5 orders of magnitude. At 
31.0 ps per cycle per cell, the solutions took about 16 hr 
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to converge on the Cray-C90. The fine grid exhibited 
slow residual and C L convergence characteristics partly 
because of the inability to increase the Courant- 
Friedrichs-Lewy (CFL) number beyond 10. The medium 
grid required 162 MW of memory and took 700 itera- 
tions to reduce the residuals by almost 2.5 orders of mag- 
nitude. The medium solution was obtained after 
approximately 4.3 hr of Cray NAS time and took 24.0 (is 
per cycle per cell. A rapid convergence was evident for 
the medium grid, and the CFL number was quickly 
ramped up to 30. The coarse grid solution was obtained 
by using the Cray-YMP because of the reduced memory 
size of 96 MW. The time per cycle on this machine was 
50.0 (is per cycle per cell, and the solution converged 
after 700 iterations, which took about 5.4 hr. As with the 
medium grid, rapid convergence was achieved for the 
coarse grid, and the CFL number was ramped up to 30.0. 
Figure 13(b) shows similar C L convergence histories for 
the medium and coarse grids. 

Figure 14 shows the convergence history for the 
adapted medium grid at a = 10° (the configuration with- 
out the air dams). The solution converged fairly quickly 
(and was ramped up to a CFL number of 30.0) after 
about 900 iterations, which took 5.2 hr on the Cray-C90. 
The convergence histories for the F-16XL with air dams 
for a= 10° and 15° are shown in figures 15 and 16, 
respectively. Compared to the adapted medium grid, the 
configuration with the air dams at a = 10° required more 
than double the number of iterations to converge. The 
poorest convergence characteristics were seen for 
a = 15°, which is shown in figure 16. After 3200 itera- 
tions, the C L was still oscillating, and a decrease of only 
one order of magnitude in the residuals was noted. The 
total run time on the Cray-C90 was 24 hr. The large 
amount of CPU time required (5 1 .7 hr if run on the Cray- 
YMP) to obtain an “almost converged” solution was 
impractical. Also the solution may have developed some 
unsteady flow-field characteristics at that point. 

Pressure and Flow-Field Analysis 

Without air dams. The fine, medium, and coarse 
grid surface pressures for the F-16XL configuration 
without the air dams at M~ 0.148 and a = 10° are com- 
pared with experimental data in figure 17 at 13 different 
fuselage stations. As stated, this comparison was done to 
determine the sensitivity of the Euler solutions (using 
USM3D) to surface grid refinement. However, the small 
benefits gained in solution accuracy do not justify the 
added memory and run time required with a finer surface 
grid. 

While the computed pressure coefficients for all the 
grids, shown in figure 17, follow the general trend of the 
experimental data, all three grids underpredicted those 


surface C P distributions influenced by vortical flow (see 
FS7.4-FS15.0). Computed C P values for all the grids 
compared better to experiment at FS16.3-FS19.7; how- 
ever, because there were no pressure ports at these sta- 
tions in the higher peak regions (except at FS19.7), the 
researcher could not determine whether the computations 
also predict suction peaks similar to the measured values. 
Compared to the medium and coarse grids, the fine grid 
solution generally yielded greater negative pressures on 
both upper and lower surfaces for all fuselage stations. 
The medium and coarse grids yielded almost identical 
pressures on the lower surface, while the fine grid pres- 
sures appeared to be offset by a negative constant. It is 
difficult to determine which grid correlated better with 
experiment on the lower surface because measured pres- 
sure data were limited; however, the relatively large 
decrease (compared to the coarse and medium grid) in 
computed lower surface C P of the fine grid appears sus- 
pect. Although pressure correlations improved slightly 
for the fine grid for the first eight fuselage stations, the 
medium and coarse grids correlated better with experi- 
ment for FS17.0-FS18.0. From this comparative analy- 
sis, it is evident that the small improvement gained from 
using the fine grid does not warrant the additional com- 
puter run time and memory associated with it. The 
medium grid was, therefore, chosen as the best surface 
grid because it captured the overall principal flow char- 
acteristics in minimal run time. 

The surface C P distributions for the adapted medium 
grid are plotted with the medium grid results in figure 18. 
The localized clustering of grid points in the regions of 
vortical flow resulted in better resolution of the off- 
surface flow field, which in turn yielded more accurate 
surface pressure distributions. The computed C P values 
obtained with the adapted grid correlated quite well with 
the experimental data, and the differences between the 
solutions for the adapted and nonadapted grids were sub- 
stantial. Both grids yielded the same lower surface pres- 
sure distributions because the point distributions were the 
same on the lower part of the wing. The computed suc- 
tion peak at FS10.3 appears to have been overpredicted, 
and a slightly lower computed suction peak is seen at 
FS13.5 and FS15.0. At the last fuselage station, the 
adapted medium grid solution showed an additional low- 
pressure region (peak at yl(bl2\ = 0.72). 

Off-body C P contours for the medium and adapted 
medium grid at FS10.0, FS18.0, and FS19.7 appear in 
figures 19, 20, and 21, respectively. Compared to the 
medium grid, the adapted grid solution yielded a more 
developed vortical flow system with larger negative pres- 
sure values in the vortex core for all fuselage stations. 
The low-pressure region at yl(bl2\ = 0.72 predicted by 
the adapted grid previously mentioned can be further 
explained by examining figure 21. The adapted vortex 
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core (fig. 21(b)), which had moved outboard of the actua- 
tor pod, affected the wing surface below it, causing the 
low suction peak seen in figure 18 at FS19.7. Because 
the primary vortex core was positioned directly over the 
actuator pod for the medium grid, only one negative suc- 
tion peak resulted, and that formed over the actuator pod 
or at y/(b/2)\ = 0.63 at FS19.7 (fig. 18). Figure 22 shows 
a close-up view of the velocity vectors in this region. The 
effect of the computed adapted vortex on the wing is evi- 
dent in figure 22(b), which shows increased circulation 
near the surface. 

The computed C P and C P t contours and crossflow 
velocity vectors obtained from the adapted grid are com- 
pared to the experimental five-hole pressure probe data 
in figures 23-31. All contour values are displayed with a 
gray scale to provide a clear depiction of the vortex size 
and shape, as well as its varying magnitude. Correspond- 
ing minimum and maximum limits at each fuselage sta- 
tion were used to scale the C P and C Pt contours. The 
computed and measured velocity vector magnitudes were 
multiplied by their appropriate nondimensional values 
for comparison purposes. Because of the size of the pres- 
sure probe, flow-field data could not be obtained very 
close to the surface. The computed C P contours are com- 
pared to experiment at FS9.0 in figure 23. The computed 
vortex-core C P value was -2.05 in the vortex core, 
located at y = 2.79 in. and z = 3.87 in., and was slightly 
inboard and closer to the wing surface than the measured 
vortex core, which had a value of -1.25 and was located 
at y = 2.82 in. and z = 3.92 in. These results appear to 
substantiate the surface pressure distributions plotted in 
figure 18 for FS8.6 and FS10.3, which show that the 
computations indeed overpredicted the experiment. 

The C P>t contours shown in figure 24 are more repre- 
sentative of the oval shape of the vortex depicted by the 
velocity vectors in figure 25. A greater total pressure loss 
is seen for the computed vortex, which had a core value 
of -1.28 compared to a -0.60 measured value. These 
numerically induced computed total pressure losses were 
largest in areas of high gradients and in regions with con- 
centrated vorticity (ref. 4). The crossflow velocity vector 
plots compared in figure 25 show similar magnitude and 
direction. The computed and measured C p contours at 
FS14.0 appear in figure 26. The computed C P for the 
vortex core located at y = 3.98 in. and z — 4.1 1 in. had a 
magnitude of -1.27, and the experimental vortex core 
positioned at y - 4. 1 1 in. and z = 4.16 in. reported a mag- 
nitude of -1.31. This trend is consistent with the C P plots 
in figure 18 at FS 13.5 and FS15.0, which show the mea- 
sured C P vortex core value to be slightly higher than the 
computed value. The vortical structures look very simi- 
lar, except at the surface where the absence of the bound- 


ary layer in the inviscid CFD solution did not affect the 
vortex shape, as seen in the experiment. The computed 
and measured vortical shapes and structures appear more 
alike in figure 27, which plots the C Pt contours. The 
computed Cpj value in the vortex center was approxi- 
mately double the measured value, which was also true at 
FS9.0. The velocity vector plots shown in figure 28 
reveal a similar pattern and trend. Both the computed and 
measured Cp values in the vortical core for station 
FS18.0 (fig. 29) were approximately -0.88. As occurred 
for the two previous fuselage stations, the predicted vor- 
tex core located at y = 5.05 in. and z = 4.22 in. was 
slightly inboard and closer to the wing surface than was 
the experimental vortex core (y = 5.25 in., z = 4.39 in.). 
The predicted and measured C P t contour plots (fig. 30) 
showed similar vortex structure, and the computed 
pressure loss was double that measured in the vortex 
core. Figure 32 shows corresponding crossflow velocity 
vectors. 

With air dams. This section presents and discusses 
the pressure distribution plots and crossflow pressure 
contours and velocity vectors for the F-16XL configura- 
tion with air dams at a = 10°, 13°, and 15°. As before, all 
grids were adapted to the vortical flow field. The compu- 
tations were compared to experiment at 10 butt-line sta- 
tions and 13 fuselage stations. 

The computational pressure distributions for a = 1 0° 
are compared with experiment in figures 32 and 33. Very 
good correlation is seen at all the butt-line stations plot- 
ted in figure 32. The extracted computational pressures at 
BL2.2 and BL2.5 are rather nonsmooth because of the 
sparse number of grid points used to define the inboard 
wing area near the fuselage. The largest discrepancy 
between the data sets appears forward on the wing at sta- 
tion BL3.8, where the computations failed to predict the 
suction peak at x/c = 0.04. This low-pressure region may 
have been caused by a secondary vortex that the Euler 
method could not capture or by leading-edge suction that 
the solution did not pick up because the grid was not fine 
enough in that area. The computational pressure results 
also correlated well with experiment at all the fuselage 
stations (fig. 33). The calculated suction peak at FS7.4 
was underpredicted and most likely the result of an insuf- 
ficient number of grid points being generated in this area. 
Compared to the configuration without the air dams in 
which only one vortex was generated from the leading 
edge, the more complicated flow with the air dam exhib- 
ited multiple vortices. The influence of the air dam on the 
flow was first seen at FS15.0, where a sharp transition 
from high to low pressure occurred at about the 0.93 
local spanwise location. A more detailed account of the 
differences between the configurations with and without 
the air dam follows. 
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The calculated particle traces for the F-16XL with 
air dams (shown in fig. 34) trace the path of each vortex, 
and are labeled for easy reference. The vortices depicted 
in the Cp contours and velocity vector plots shown in 
figures 35-37 are labeled likewise. Figure 35 shows the 
vortex generated from the apex of the leading edge in a 
crossflow plane at FS10.0. This vortex, labeled A, propa- 
gated along the inboard portion of the wing. A smaller, 
clockwise-rotating vortex, B, emanated from the sharp 
leading edge of the air dam and traced an outboard path 
(fig. 36). Another vortex, C, which rotated in the counter- 
clockwise direction and appeared outboard of vortex B, 
formed as a result of a saddle-point or off-surface separa- 
tion in the flow. The effect of this vortical system on 
the surface is evident in the pressure distribution plot 
at FS16.3 (fig. 33). Vortex A had a suction peak at 
y/(b/2)i = 0.72, and vortices B and C had C P values of 
approximately -0.60 and -0.70 at yl(bl2\ = 0.84 and 
0.89, respectively. As the vortical system propagated 
downstream, it continued to expand, and a newly formed 
vortex D appeared at FS20.2 (fig. 37). This vortex must 
have been very weak because there is no intertwining of 
the particle traces in figure 34. Figure 38 shows a close- 
up view of the velocity vectors in the missile area at 
FS20.2. One sees a vortex on three of the missile fins and 
one in the middle of the upper two fins. 

The computed versus experimental surface pressure 
coefficient comparisons at a=13° appear in figures 39 
and 40. The computational pressures match well at this 
angle of attack, and the grid, which was adapted for the 
cl - 10° solution, appears sufficiently clustered to accu- 
rately resolve the flow at a = 13°. Figures 41 — 43 show 
the computed pressure contours and velocity vectors in 
three crossflow planes. The overall flow features and 
vortical footprints were similar to the a= 10° solution; 
therefore, the off-surface particle traces are not pre- 
sented. Examination of the labeled vortices depicted in 
figures 4 1 — 43 reveals the same vortical system identified 
for oc = 10°. At FS10.0, one sees a single vortex, A, on 
the wing that is similar (although greater in magnitude) 
to the one that developed for the a = 10° solution. Like- 
wise, at FS16.3, shown in figure 42, multiple vortices 
developed as a result of the flow impinging on the air 
dam, and vortex A appears to have split in two. However, 
the off-surface particle traces revealed only a single vor- 
tex in that region, and as shown in figure 43, vortex A 
did indeed remain intact. In contrast to the a = 10° solu- 
tion, vortex C appears to be merging with the others in 
the flow field, and vortex D is more developed. The 
velocity vectors plotted around the missile in a crossflow 
plane are shown in figure 44. A small vortex, which did 
not exist for the a = 10° solution, formed on the lower 
left fin. 


The computed pressure distributions at all fuselage 
and butt-line stations correlated well with experiment for 
a = 15°, as shown in figures 45 and 46. Compared with 
the a = 10° and a = 13° solutions, which showed a 
smooth distribution of pressures on the upper surface at 
the last three butt-line stations, the a =15° solution 
yielded a more uneven distribution of pressures (BL5.9- 
BL6.9 in fig. 45). This uneven pressure distribution may 
have been caused by the increase in spanwise flow on the 
outboard portion of the wing. In general, computational 
and experimental pressure distribution trends remained 
similar for all three angle-of-attack solutions. The com- 
puted off-surface particle traces are shown in figure 47, 
and the computed pressure contours and velocity vectors 
are presented in three crossflow planes in figures 48-51. 
Unlike the previous two solutions in which vortex B 
remained intact near the trailing edge, vortex B merged 
with vortex A at FS20.3 (fig. 50). Also vortex C occurred 
more outboard on the wing as a result of the increased 
spanwise flow. In figure 50, the increased circulation of 
the vortex A,B and vortex C appears to prevent the for- 
mation of the smaller vortices evident in the lower angle- 
of-attack solutions. Flow in the vicinity of the wingtip 
missile fins shown in figure 5 1 closely resembled that of 
the a = 13° solution (fig. 39) only more well-defined. 

Figure 52 shows the effect of adding the air dam on 
the overall surface Cp of the F-16XL for oc = 10°. The Cp 
pattern remained basically the same for both configura- 
tions up to about FS13.5, and similar computational and 
experimental C P data were noted at this station (fig. 53). 
In contrast to the surface flow without the air dam, the 
vortical flow of the air-dam configuration split and was 
directed to each side of the air dam, as shown in 
figure 52(b). The flow that traveled outboard accelerated 
spanwise toward the wingtip, which resulted in greater 
negative C P values on the outboard wing compared to 
those of the no-air-dam configuration. Higher suction 
peak values were noted for both experiment and compu- 
tations (outboard of the air dam) for FS16.3-FS19.7, 
shown in figure 53. The experimental peak pressure val- 
ues at FS13.5 and FS15.0 were slightly higher for the 
model without the air dam (fig. 53) although the compu- 
tations show the opposite. This discrepancy probably 
was caused by the grid. Because more points were 
required to resolve the air dam, the higher grid density in 
this area may have resolved the flow better. 

Concluding Remarks 

The subsonic (M = 0.148) flow field about the basic 
F-16XL cruise configuration, with and without air 
dams, was modeled with an inviscid unstructured grid 
technique, namely VGRID and USM3D. In order to 
assess the grid dependency of the solution, a grid study 
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was performed in which coarse, medium, and fine grids 
were generated for the model without air dams. No sig- 
nificant improvements in solution accuracy were seen 
(through comparison of computations to experiment) 
when the number of grid points was globally increased. 
In fact, all the USM3D solutions that used these grids 
exhibited poor correlation with experiment when com- 
pared to the adapted grid. The finest grid, which had over 
1000000 points, exhibited very slow convergence char- 
acteristics and took about 16.0 hr on the Cray-C90 to 
converge. The grid distribution of the medium grid was 
chosen as the most efficient surface representation in 
terms of memory, run time, and convergence qualities. 
The researcher used the medium grid as a starting point, 
and the off-surface grid of the basic configuration with 
air dams was locally adapted to the vortical flow at 
a= 10°, 13°, and 15°. Increases in the grid-point density 
in the vortical regions resulted in a more resolved flow 
and a marked improvement in the comparison between 
computational and experimental pressure distributions 
when compared to the nonadapted grid. An integrated 
package that couples the adaptive grid process with the 
grid generator and flow solver would be an invaluable 
tool for the designer and researcher. Not only would the 
grid generation time and effort be reduced but also any 
guess work on the part of the investigator in terms of 
locating the vortex core for adaption purposes would be 
eliminated. 

The computed off-surface static and total pressure 
coefficient contours and velocity vectors compared well 
with the BART five-hole probe flow-field data (F-16XL 
without air dams) and provided a more comprehensive 
description of the flow physics when used in conjunction 
with the surface pressure distributions. For all crossflow 
planes examined, the computed vortex-core static- 
pressure-coefficient value was located slightly inboard 
and closer to the wing surface than was the measured 
value. In contrast to the F-16XL model without air dams, 
in which only one vortex was generated from the leading 
edge, the flow with the air dams was more complicated 
in that multiple vortices were produced. Comparison of 
computed results with experimental data shows that most 
of the primary features of the flow were correctly simu- 
lated. Differences between the predicted and measured 


results were attributed mainly to viscous effects, such as 
secondary separation. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
August 26, 1996 
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(a) Coarse grid. 




(c) Fine grid. 

Figure 6. Surface grids of varying densities for F-16XL without air dams. 
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(a) Coarse grid. 


(b) Medium grid. 




(c) Fine grid. 

Figure 8. Cross-sectional grid planes at FS18.0 for F-16XL without air dams. 
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(a) Medium grid. 



(b) Adapted medium grid. 

Figure 9. Adapted and nonadapted surface grids for F-16XL without air dams. 







(a) Planform view of upper surface mesh. 



(b) Close-up view of air-dam geometry. 


Figure 12. Adapted surface mesh for F-16XL with air-dams, a = 10°. 





Log (residual) Log (residual) 



(a) Residual convergence. (b) C L convergence. 

Figure 14. Convergence characteristics of adapted medium grid for F-16XL without air dams, a = 10°, M = 0.148. 



21 










o BART upper surface 

• BART lower surface 

CFD fine 

CFD medium 

CFD coarse 




Figure 17. Computed spanwise C P distributions from grid study compared to experiment for F-16XL without air dams, 
oc = 10°, M = 0.148. 
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Figure 17. Concluded. 
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Figure 18. Computed spanwise C P distributions from adapted and nonadapted grids compared to experiment for 
F-16XL without air dams, a — 10 °,M = 0.148. 
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Figure 18. Concluded. 









(a) Medium grid. 



(b) Adapted medium grid. 

Figure 22. Computed velocity vectors for adapted and nonadapted grids at FS19.7 for F-16XL without air dams, 
a = 10°, M = 0.148. 
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(a) Computed C P contours. (b) Experimental C P contours. 

Figure 23. Computed and measured C P contours at FS9.0 for F-16XL without air dams, a = 10°, M 


(a) Computed C P contours. (b) Experimental C P contours. 

Figure 24. Computed and measured C P t contours at FS9.0 for F-16XL without air dams, a = 10°, M 
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(a) Computed velocity vectors. 
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(a) Computed C P contours. (b) Experimental Cp contours. 

Figure 26. Computed and measured C P contours at FS14.0 for F-16XL without air dams, a = 10°, M 


(a) Computed C P t contours. (b) Experimental C P t contours. 

Figure 27. Computed and measured C P t contours at FI 4.0 for F-16XL without air dams, a = 10°, M 









(b) Experimental velocity vectors. 


Figure 28. Computed and measured velocity vectors at FI 4.0 for F-16XL without air dams, a 






(a) Computed C P contours. (b) Experimental C P contours. 

Figure 29. Computed and measured C P contours at FS18.0 for F-16XL without air dams, a = 10°, M — 0.148. 



(a) Computed C P t contours. (b) Experimental C P t contours. 

Figure 30. Computed and measured C P t contours at FS 18.0 for F-16XL without air dams, a = 10°, M = 0.148. 
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(a) Computed C P contours. 



(b) Close-up view of computed velocity vectors. 

Figure 37. Computed Cp contours and velocity vectors at FS20.2 for F-16XL with air dams, a = 10°, M = 0.148. 
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Figure 38. Close-up view of missile region depicting computed velocity vectors at FS20.2 for F-16XL with air dams, 
10°, M = 0.148. 
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(a) Computed Cp contours. 



(b) Computed velocity vectors. 


Figure 41. Computed C P contours and velocity vectors at FS 10.0 for F-16XL with air dams, a = 13°, M = 0.148. 
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(b) Close-up view of computed velocity vectors. 

Figure 42. Computed C P contours and velocity vectors at FS 16.3 for F-16XL with air dams, a - 13°, M = 0.148. 
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(a) Computed C P contours. 



(b) Close-up view of computed velocity vectors. 

Figure 43. Computed C P contours and velocity vectors at FS20.2 for F-16XL with air dams, a = 13°, M = 0.148. 
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Figure 44. Close-up view of missile region depicting computed velocity vectors at FS20.2 for F-16XL with air dams, 
a = 13°, M = 0.148. 
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(a) Computed C P contours. 



(b) Computed velocity vectors. 

Figure 48. Computed C P contours and velocity vectors at FS10.0 for F-16XL with air dams, a 
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Figure 53. Span wise distributions of experimental and computational C P for F-16XL with and without air dams, 
a = 10°, M = 0.148. 
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